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Knowing  the  Past  is  Good,  but  Knowing  the 
Future  is  even  Better:  Predictive  Models! 
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First-Principles  Models  (physiology  based) 


Data-Driven  Models  (derived  from  data) 
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Ordinary  Differential  Equation 
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Model  Requirements  for  Practical  Use 


•  Highly  accurate  (for  a  reasonable  horizon) 

^  Individual-specific  models 

•  Minimum  manual  tuning 

^  Adaptive  first-principles  models 
^  “Universal ’’data-driven  models 

-  need  to  measure  ^‘something”  from  the  individual  - 


Physiologic  Variables  (or  States) 
We  Wish  to  Predict 


1 .  Performance  impairment  due  to  totai  sieep  ioss 

2.  Body  core  temperature  (minimize  heat  injuries) 

-  Armed  forces  (2003-2005)^:  3617  heat  exhaustion  injuries 

784  heat  stroke  injuries 

3.  Giucose  concentrations  (diabetic  patients) 

-  $1  in  every  $10  heaith  care  doiiars  is  attributed  to  diabetes* 


■•■chppm-www.apgea.army.mil/heat/  *Hogan  et  al.,  Diabetes  Care,  26,  917  (2003) 


PVT  lapses  (greater  impairment  -^) 


Performance  Impairment  Prediction 

Problem:  Predict  performance  impairment  due  to  total  sleep 
loss  (82-hour  study  of  total  sleep  deprivation)* 

Measure  of  Performance:  Lapses  in  a  reaction-time  test 
(Psychomotor  Vigilance  Task  -  PVT)  every  two  hours 


*Wesensten  et  al,  J.  Sleep  Res,  14,  255  (2005) 

PVT  lapses:  #  of  times  reaction  time  >  500  msec  over  a  10-min  test  session 


Data  from  Tom  Balkin’s  group  (Walter  Reed) 


Two-Process  Model  of  Sleep  Regulation 

-  “state-of-the-art”  first-principles  model  - 


Process  S 

(sleep/wake  history) 

S(t)  =  1  -  exp(-Ar  /  r.)(l  -  S(t  - 1)) 
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'Borbely,  Human  Neurobiol.,  1,  195  (1982) 


Predicting  Performance  Using  the 

Two-Process  Model 


Performance  P(t): 

P(t)  =  S(t)  +  C(t) 

For  total  sleep  deprivation: 

5 


P(t)  =  a  -  05(0) 


Population-Avg.  Solution:  Model  parameters  (r^,  S(0),  a,  are  fixed 


PVT  lapses  (greater  impairment 


Population-Avg.  Model  Not  Optimal 


Solution:  Individual-specific  adaptive  models 

-  Models  parameters  (r^,  S(0),  a,  (5,  )  are  automatically 

adjusted,  for  each  individual,  after  each  PVT  observation 


First-Principles,  Two-Process  Model 

-  Adaptive,  Individual-Specific  Model  - 

Sleep-Deprived 
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Adaptive  Individual-Specific  Models* 


Performance  (P): 
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*  Rajaraman  et  al.,  J  AppI  Physiol.,  104,  459  (2008) 


PVT  lapses  (greater  impairment 


Predicting  a  Vulnerable  Subject 
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PVT  lapses  (greater  impairment 
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PVT  lapses  (greater  impairment 


Individualized  vs. 

Predictions 


2-hour  ahead  prediction 


'PVT  lapses 

'Population  average  (RMSE  =  18.9) 
'Individualized  (RMSE=  14.3) 


Vulnerable 
(subject  #44) 

J ^ ^ ^ ^ ■  ■  ■ 
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Rajaraman  et  al.,  J  AppI  Physiol.,  104,  459  (2008) 


PVT  lapses  (greater  impairment 


Predicting  a  Resilient  Subject 


PVT  lapses  (greater  impairment 
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PVT  lapses  (greater  impairment 


Individualized  vs. 

Predictions 


2-hour  ahead  prediction 
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Rajaraman  et  al.,  J  AppI  Physiol.,  104,  459  (2008) 


Autoregressive  (AR),  Data-Driven  Model 


y{t)  =  b^y(t  -  1J  +  b2y{t  -  2)  +  ...  +  b^oy{t-  30) 
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y{t)\  measurement/prediction  at  time  t 


b{.  model  coefficients  (unknown) 


30:  number  of  previous  measurements 


Found  to  be  individual 
independent!!! 


Autoregressive  (AR),  Data-Driven  Model 

-  “Universal,  ”  Individual-Independent  Model  - 


Continuous  Glucose 
Monitor 
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Data-Driven, 
Universal 
Model  (fixed 
coefficients  b:) 


Core  Temperature, 


Body  Core  Temperature  Prediction 


Problem:  Predict  core  temperature  20-minutes  ahead’ 
Measurement:  Temperature  pill  (1 -minute  data) 
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Joint  Readiness  Training 
Center,  Ft.  Polk,  LA,  Aug/01 

•  West  Point  Cadet  and 
Soldier  from  509th  Infantry 
Regiment 

•  Soldier  was  very  fit  and 
heat  acclimated 

•  Cadet  was  not  as  fit,  not 
heat  acclimated  and  was 
carrying  a  heavier  load 

•  Soldier’s  model  used  to 
predict  Cadet’s  temperature 
and  vice  versa 


Data  from  Reed  Hoyt’s  group,  U.S.  Army  Research  Institute  of  Environmental  Medicine  (USARIEM) 


Core  Temperature, 


Cadet’s  Model  Used  to  Predict  Soldier’s 
Core  Temperature  20-min  Ahead 
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Soldier’s  Model  Used  to  Predict  Cadet’s 
Core  Temperature  20-min  Ahead 
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simulated  real-time  predictions 


Gribok,  Duller,  and  Reifman,  IEEE  Trans.  Biomed.  Eng. ,  in  press 
Gribok,  Duller,  Hoyt,  and  Reifman,  under  review 


Glucose  (mg/dl) 


Glucose  Prediction  for  Type  1  &  2  Diabetes 

-  three  studies  using  distinct  continuous  glucose  monitoring  (CGM)  devices  - 
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*Data  provided  by  Ken  Ward  (iSense),  Robert  Vigersky  (DexCom),  Direcnet  (Guardian  RT) 


Glucose  Prediction  for  iSense  Subject  #7 


Raw  Data 


Reifman  et  al.,  Diabetes  Science  &  Technology,  1,  478  (2007) 
Gani  et  al.,  under  review 


Models  from  3  Different  Subjects  Used  for 
Glucose  Prediction  of  iSense  Subject  #7 
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Universal  Model 

-  predictions  across  different  subjects,  devices,  types  of  diabetes  - 


30-nnin  Prediction  Horizon  (models  based  on  different  subjects/devices) 


Testing 

iSense  Model  (Avg.  8) 

Guardian  RT  Model  (Avg.  18) 
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Gani  et  al.,  in  preparation 


What’s  Next?  Field  Testing  of  Universal  Model 


Real-Time  Body  Core 
Temperature  Prediction 


y{t)  =  b^y{t  - 1)  +  b2y(t  ■  2)  +  ...  +  bs^yit  -  30) 


Joint  effort  with  Reed  Hoyt’s  group  (USARIEM) 


“All  models  are  wrong, 
some  are  useful.  ” 


George  Box 


BkJtectinology  HPC 
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